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ABSTRACT 

The gravitational evolution of the genus and other statistics of isodensity contours 
of the density field is derived analytically in a weakly nonlinear regime using second- 
order perturbation theory. The effect of final smoothing in perturbation theory on 
the statistics of isodensity contours is also evaluated. The resulting analytic expres- 
sion for the genus is compared with A^-body numerical simulations and exhibits a 
good agreement. 



1. Introduction 

There are many statistical tools to investigate the three-dimensional pattern of the 
galaxy distribution including two-point correlation function, power spectrum, cluster 
correlations, higher-order correlations, probability distribution function, etc. Among 
others, there is a class of statistics using a smoothed density field which cut the noisy 
property of galaxy distribution. 

For example, the topology^ of isodensity contours of those smoothed field is the 
popular one. The genus G, which is defined by —1/2 times the Euler characteristics of 
two-dimensional surfaces, divided by the total volume, can be a quantitative measure 
of the topology. The genus is a function of smoothing scales and the density threshold. 
The genus as a function of density threshold for a fixed smoothing scale is called the 
genus curve and is analyzed both in numerical simulations and in redshift surveys of 
galaxies by many people^. 

There are other statistics of isodensity contours, which include the 2D genus in 
two-dimensional slices of density field"^, the area of isodensity contours'*, A3, the 
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length of isodensity contours in two-dimensional slices'*, and the level crossing 
statistics"*, A^i. 

The analytical expression of these statistics are known for Gaussian random den- 
sity field"^''*'^. They are 

(1) 
(2) 

2Nr{v) = -iV2(^) = iV3(^) = - f e-^'^\ (3) 

TT IT \ 6 J 

where p is the difference between density threshold and mean density in units of 
standard deviation of density and 

with P{k) being the power spectrum of the density fluctuation. 

Previous analyses mainly compared the observed genus etc. with the random 
Gaussian prediction (f )-(3). With sufficiently large smoothing scales, this comparison 
could tell us if initial density fluctuation is random Gaussian or not. With the finite 
smoothing scale of cosmological interest, the effect of nonlinear gravitational evolution 
would be substantial. This nonlinear effect has been explored only by using A^-body 
numerical simulations so far. In this paper, I present the analytic expression of 
statistics of isodensity contours extending the result for genus in the weakly nonlinear 
regime using a second-order perturbation theory^. 



2. Isodensity Statistics for Quasi-Gaussian Random Field 

First, I introduce the seven quantities for a non-Gaussian random field 6[x^y^z) 
with zero mean as (A^) = (7~^{6^d6 / dx^ dS / dy ^ 86 / dz^ / dx'^ ^ / di/ ^ / dxdy) 

where a = \J (S"^) is rms of the field and the field is defined in Cartesian coordinates 
x^y^z. The statistics of isodensity contour (7,(72,^1,^2,^3 of constant surfaces 
6 = per is given by"*'^ 

G{p) = {Sd{Ai - p)SD{A2)SD{As)\A,\{AsAe - A^)) , (5) 

G2{'y) = -\ {8{A, - v)8d{A,)\AM^) , (6) 

2Nr{v) = -N2{v) = N^{v) = 2 {8{Ar - v) \A^\) , (7) 

TT 

where 8]^ is a Dirac's delta-function. This expression is valid for general homogeneous 
and isotropic non-Gaussian random fields. 
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In the following, I assume that (A^j • • • ~ 0(cr^-2). This relation IS a very 

definition of "weak non-Gaussianity" in this paper and is a result of perturbation 
theory''. Using the multi-dimensional version of the Edgeworth expansion^ and similar 
calculation developed by Matsubara^, the above expression is expanded to the first 
order in a as follows: 



G{i 



1 



(27r] 



G2{v) : 
2iVi(//) 



1 (fe^) -.V2 

(27r)3/2 3 



(8) 
,(9) 



TT 
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where Hn{v) = ( — )"e'' ' {dl dv)'^e ^ ' are Hermite polynomials, and I have defined 
three quantities, 

1 



T 

U 



1 



2(F)(7 
3 



11^ 



4(F)V' 



-(V<5- V<5V'<5). 



The quantity S is usually called "skewness" . The first term in square brackets of 
Eqs. (8)-(10) corresponds to Gaussian contribution and the other terms correspond 
to non-Gaussian contribution. 



3. Gravitational Evolution of the Genus Curve in Second Order Pertur- 
bation Theory 

Gravitational nonlinear evolution give rise to T, U even from the initial Gaus- 
sian random density fiuctuation which has vanishing S , T, U . I use second order 
perturbation theory of the non-relativistic coUisionless self-gravitating system in the 
fiuid limit^ to compute T, U to lowest order in a in Einstein-de Sitter universe. 
The result including the effect of smoothing evolved field by Gaussian filter is given 
by Matsubara^. In fact, the observable curve is obtained by smoothing of density 

^The usual one-dimensional Edgeworth expansion was recently applied to gravitational instability 
theory® 
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unsmoothed 


n = 1 


n = 


n = —1 


n = -2 


n = —3 


CDM 


HDM 


s 


4.857 


3.029 


3.144 


3.468 


4.022 


4.857 


3.443 


4.018 


T 


3.905 


2.020 


2.096 


2.312 


2.681 


3.238 


2.335 


2.970 


U 


1.543 


1.431 


1.292 


1.227 


1.222 


1.272 


1.223 


1.306 



Table 1. The numerical values of S*, T, U for unsmoothed perturbation theory and smoothed 
perturbation theory for power-law spectra, n = 1 to — 3, CDM and HDM models. 



CD 




Fig. 1. Comparison of the analytic genus curve (a solid line) with that from an A^-body simulation 
(symbols) for the Poisson model at (7 = 0.2. A dashed line indicate the Gaussian prediction. 



fluctuation at the final stages. Numerical values of S*, T, U for power-law initial spec- 
tra, CDM and HDM models are summarized in Table 1. In CDM and HDM models, 
I set Oo = 1, Hq = 50 km s~^ Mpc~^ and smoothing length = 10 Mpc. The values of 
skewness S in this table was first obtained by Juszkiewicz et al.^° 

4. Comparisons with Numerical Simulations 

Suto and I found that the analytic expression of the genus curve agrees well with 
the numerical A^-body simulations for power-law initial fiuctuations P{k) oc and 
CDM models^^. In Fig. 1 is plotted the normalized genus G{p) / G{0) as an example 
of the comparison. The initial fiuctuation is the Poisson model P{k) = const. The 
rms of the fiuctuation at the stage of comparison is cr = 0.2. 

5. Discussion 

The prominent feature of the results, Eqs. (8)-(10) is that weakly non-Gaussian 
correction introduces asymmetry to the symmetric or anti-symmetric curves. The 
pattern of the asymmetry is dependent on initial power spectra through smoothing 
effect of S*, r, U . Thus, in principle, accurate observations on the statistics presented 
here can restrict the properties of initial fiuctuation, such as Gaussianity, the shape 



4 



of the spectrum, by the amphtude and the pattern of asymmetry of the curve. The 
projects as Sloan Digital Sky Survey (SDSS) will enable us to have a large amount of 
redshift data in near future and the analysis indicated in this paper will be important. 

6. Acknowledgements 

I am grateful to thank Y. Suto for useful discussions. I acknowledge the support 
of JSPS Fellowship. This research was supported in part by the Grants-in-Aid for Sci- 
entific research from Ministry of Education, Science and Culture of Japan (No. 0042). 

7. References 



1. J. R. Gott, A. L. Melott and M. Dickinson, Astrophys. J. 306 (1986) 341. 

2. See, e.g., M. S. Vogeley, C. Park, M. J. Geller, J. P. Huchra and J. R. Gott, 
Astrophys. J. 420 (1994) 525; and references therein. 

3. A. L. Melott, A. P. Cohen, A. J. S. Hamilton, J. R. Gott and D. H. Weinberg, 
Astrophys. J. 345 (1989) 618. 

4. B. S. Ryden, Astrophys. J 333 (1988) L41; 

B. S. Ryden, A. L. Melott, D. A. Craig, J. R. Gott, D. H. Weinberg, R. J. Scher- 
rer, S. P. Bhavsar and J. M. Miller, Astrophys. J. 340 (1989) 647. 

5. R. J. Adler, The Geometry of Random Fields (Wiley, Chichester, 1981); 

A. G. Doroshkevich, Astrophysics 6 (1970) 320 (transl. from Astrofizika 6 
581); 

J. M. Bardeen, J. R. Bond, N. Kaiser and A. S. Szalay, Astrophys. J. 304 
(1986) 15; 

A. J. S. Hamilton, J. R. Gott and D. H. Weinberg, Astrophys. J. 309 (1986) 

1; 

J. R. Bond and G. Efstathiou, Mon. Not. R. Astron. Soc. 226 (1987) 655. 

6. T. Matsubara, Astrophys. J. L43 (1994) L43. 

7. J. N. Fry, Astrophys. J. 279 (1984) 499; 

M. H. Goroff, B. Grinstein, S.-J. Rey and M. B. Wise, Astrophys. J. 311 (1986) 
6; 

F. Bernardeau, Astrophys. J. 392 (1992) 1. 

8. R. Juszkiewicz, D. H. Weinberg, P. Amsterdamski, M. Chodorowski and 
F. R. Bouchet, Astrophys. J. (1994) in press; 

F. Bernardeau and L. Kofman, IFA-94-19 (1994) preprint. 

9. P. J. F. Peebles, The Large-Scale Structure of the Universe (Princeton, Prince- 
ton University Press, 1980) §18. 

10. R. Juszkiewicz, F. R. Bouchet and S. Colombi, Astrophys. J. 412 (1993) L9. 

11. Y. Suto and T. Matsubara, in preparation 



5 



